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ABSTRACT 

We calculate the evolution of linear density and temperature perturbations in a 
universe with dark matter, baryons, and radiation, from cosmic recombination until 
the epoch of the first galaxies. In addition to gravity, the perturbations are affected 
by electron scattering with the radiation, by radiation pressure, and by gas pressure. 
We include the effect of spatial fluctuations in the baryonic sound speed and show 
that they induce a > 10% change in the baryonic density power spectrum on small 
scales, and a larger change on all scales in the power spectrum of gas temperature 
fluctuations. A precise calculation of the growth of linear perturbations is essential 
since they provide the initial conditions for the formation of galaxies and they can 
also be probed directly via cosmological 21cm fluctuations. We also show that in 
general the thermal history of the cosmic gas can be measured from 21cm fluctuations 
using a small-scale anisotropic cutoff due to the thermal width of the 21cm line. 
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1 INTRODUCTION 

Observations of the cosmic microwave background (CMB) 
show that the universe at cosmic recombination (redshift 
a ~ 10 3 ) was remarkably uniform apart from spatial fluctua¬ 
tions in the energy density and in the gravitat ional potential 
of roughly one part in 10 s (e.g., iBennett et all Il99ft) . The 
primordial inhomogeneities in the density distribution grew 
over time and eventually led to the formation of galaxies as 
well as galaxy clusters and large-scale structure. 

Different physi cal processes contributed to the pertur- 
bation growth (e.g., |Peeblesll980Ld993 : Ma fc Bertschingerl 
li 99 Fill . In the absence of other influences, gravitational forces 
due to density perturbations imprinted by inflation would 
have driven parallel perturbation growth in the dark mat¬ 
ter, baryons and photons. However, since the photon sound 
speed is of order the speed of light, the radiation pressure 
produced sound waves on a scale of order the horizon and 
suppressed sub-horizon perturbations in the photon den¬ 
sity. The baryonic pressure similarly suppressed perturba¬ 
tions in the gas below the much smaller baryonic Jeans 
scale. Since the formation of hydrogen at recombination had 
decoupled the cosmic gas from its mechanical drag on the 
CMB, the baryons subsequently began to fall into the pre¬ 
existing gravitational potential wells of the dark matter. 

* E-mail: smadar@wise.tau.ac.il (SN); barkana@wise.tau.ac.il 
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Spatial fluctuations developed in the gas temperature 
as well as in the gas density. Both the baryons and the 
dark matter were affected on small scales by the temper¬ 
ature fluctuations through the gas pressure. Compton heat¬ 
ing due to scattering of the residual free electrons (con¬ 
stituting a fraction ~ 10 -4 ) with the CMB photons re¬ 
mained effective, keeping the gas temperature fluctuations 
tied to the photon temperature fluctuations, even for a time 
after recom bi natio n. In prior analyses (e.g., iPeebles 1980 : 
iMa &Be^chinge 3ll995|) a nd in the standard CMBFAST 
code dSebaj^^jaMarriagalllOOflll a spatially uniform speed 
of sound was assumed for the gas at each redshift. This 
assumption meant that the gas temperature distribution 
at any given time was assumed to be simply proportional 
to its density distribution. While this assumption can be 
accurately invoked before recombination, when the gas is 
thermally and mechanically coupled to the radiation, and 
in the absence of heating would be accurate after recom¬ 
bination (i.e., for an adiabatic gas), this assumption is in 
fact inaccurate during the post-recombination era due to 
the continual Compton heating. A full calculation of the 
temperature fluctuations is essential for obtaining accurate 
results for the small-scale density power spectra of both the 
baryons and (to a lesse r deg ree) the dark matter (se e also 
lYamamoto et ahl lll99it) and lYamamoto et alJ 1^998), who 
considered the effect of a varying sound speed on density 
perturbations). Note that whether we discuss perturbations 
driven by gravity or sound waves driven by pressure gradi- 
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ents, we in every case refer to Sp/Sp as the square of the 
sound speed of the fluid, where 5p and 5p are the pressure 
and density perturbations, respectively. 

The primordial inhomogeneities in the cosmic gas in¬ 
duced variations in the optical depth for absorption of the 
CMB at the redshifted 21cm wavelength of neutral hydro¬ 
gen. Therefore, the hyperfine spin flip transition of neutral 
hydrogen is potentially the most promising tracer of the cos¬ 
mic gas in the era before the first galaxies. Future observa¬ 
tions of the redshifted 21cm power spectrum as a function 
of wavelength and angular direction should provide a three- 
dimensional map of the distribution of neutral h ydrogen 
(’Hogan fc Reedll979t iMadan. Meiksin fcReejll997th where 
redshift supplies the line-of-sight distance. At different red- 
shifts, different signal characteristics are expected, accord¬ 
ing to the evolution of density fluctuations and the thermal 
hist ory of the cosmic gas. 

iRa.rkana. fa Loeb (I2005all explored the consequences of 
the varying baryonic sound speed for the evolution of pertur¬ 
bations on sub-horizon large-scale structure scales (wavevec- 
tors k = 0.01-40 Mpc -1 ) in the redshift range z = 20-150, 
where the analysis involves only gravity; they showed that 
observations of 21cm fluctuations at these redshifts can mea¬ 
sure five easily-separated fluctuation modes. In this paper 
we extend the redshift range back to cosmic recombination, 
we consider the largest scales (§2) where general relativity 
and photon fluctuations play a role, and we also study the 
smallest scales (§3) where gas pressure has an essential ef¬ 
fect. This allows us to study the evolution of temperature 
fluctuations on all scales, and to produce the first accurate 
predictions at k > 100 Mpc -1 of the baryonic density power 
spectrum and (in §4 & 5) of 21cm fluctuations. We also cal¬ 
culate (§5) a general feature of the 21cm fluctuation power 
spectrum, i.e., smoothing by the thermal width of the 21cm 
absorption line. We give our conclusions in §6. 


2 GROWTH OF LARGE SCALE DENSITY 

PERTURBATIONS 

Up to recombination the baryons and the photons were 
tightly coupled, and afterwards the baryons were signifi¬ 
cantly Compton heated down to z < 200 and retained 
some memory of this heating at even lower redshifts. Thus, 
calculations of the growth of perturbations in the baryon 
density and temperature must take into consideration the 
evolution of the photon density and temperature fluctua- 
tions. The perturbations in the photons themselves (e.g., 
iMa fc BertschingerHl995 f are not significantly affected by 
the baryon temperature since the photon pressure dominates 
strongly over the baryonic pressure. 

If we define the photon temperature T-, at each point 
as the temperature averaged over all photon directions, then 
dimensionless fluctuations in this temperature (<5t 7 ) are re¬ 
lated to photon density fluctuations (S y ) by 

s ~t = 4<5t 7 ■ (1) 

We can derive an equation for the evolution of fluctuations 
5t in the gas temperature T from the first law of ther- 
modyn amics [following the derivation by iBarkana fc Loebl 
( 2005(1) 1: 


dQ = ^ksdT — k B Td\ogph , (2) 

where dQ is the heating rate per particle and pb is the baryon 
density. In the post-recombination era before the formation 
of galaxies, the only external heating arises from Thomson 
scattering of the remaining free electrons with CMB pho¬ 
tons, resulting in a heating rate per particle 

^ = 4^2— ksiTj - T)p lXe {t) , (3) 

at m e 

where <tt is the Thomson scattering cross section, p 7 is the 
photon energy density, and x e (t) is the electron fraction out 
of the total number density of gas particles at time t. After 
cosmic recombination, x e (t) changes due to the slow recom¬ 
bination rate of the residual ions: 

^ = -a B (T)x 2 e {l + y)n H , (4) 

where tlh is the mean number density of hydrogen at time 
t, y = 0.079 is the helium to hydrogen number density ra¬ 
tio, and a B (T) is the case-B recombination coefficient of 
hydrogen. We assume that the residual electron fraction is 
uniform. Fluctuations in x e are expected to be very small, 
since the gas in every region starts out fully ionized and 
the final x e depends on the total accumulated number of re¬ 
combinations; most recombinations occur early during cos¬ 
mic recombination, and the recombination rate slows down 
greatly once most of the gas has recombined, making x e in¬ 
sensitive at that point to the gas density and temperature. 
E.g., since the cosmic mean x e declines by only a factor 
of 1.3 between z = 100 and 20, we estimate that in this 
redshift range, fluctuations in x e are smaller than density 
fluctuations by a factor of ~ 20. 

Combining equations and 0 we obtain 


dT = 2^_, dlog pb + x e (t) p 7 _ T ) a -4 


dt 3 dt 
where 

1 _ 8 _ 0 or c 
7 “ 3 p7 m e 




= 8.55 x 10 -13 yr -1 


(5) 

( 6 ) 


and a is the scale factor. Thus the evolution of the cosmic 
mean gas temperature is 

^ = _2RT+^(f 7 -T)a- 4 , (7) 


where T 7 = [2.725 K]/a is the mean CMB temperature, and 
the first-order equation for the perturbation is 


dSr 

dt 


2dtfb ®e(t) _ 4 f i 

3 dt t-y 1 V T 




.( 8 ) 


The first term on the right-hand-side of each of these two 
equations accounts for adiabatic expansion of the gas, and 
the remaining terms capture the effect of the thermal ex¬ 
change with the CMB. 

We have numerically calculated the evolution of 
the perturbations by mo difying the CMBFAST code 
iSeliak fc"Zaldarriagajl996il according to these equations for 
the temperature perturbations (along with other modifica¬ 
tions on small scales that are described in the next section). 
The initial conditions for the temperature at z > 10000 are 
St = because of the tight thermal coupling between the 
gas and the CMB. At lower redshifts, the mechanical cou¬ 
pling weakens but the thermal coupling remains strong, i.e., 
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Figure 1. Power spectra of density and temperature fluctu¬ 
ations vs. comoving wavenumber, at redshifts 1200, 800, 400, 
and 200. We consider fluctuations in the CDM density (solid 
curves), baryon density (dotted curves), baryon temperature 
(short-dashed curves), and photon temperature (long-dashed 
curves). 


the coefficient of the coupling term in equation 0 remains 
very large, and it is numerically highly inefficient to solve 
this equation directly. We instead develop a simple approx¬ 
imation that may be used during this era of thermal tight 
coupling. First we write the expression for d(5r — ST-,)/dt 
using equation 0 along with the expression for dSr^/dt 
as computed by CMBFAST. As long as the thermal cou¬ 
pling is effective, the difference St — St 7 is very small, i.e., 
\St~ St. 7 | <S |<5b |, so that the various terms in the expression 
for 5(5t — St~ )/dt must cancel each other nearly completely. 
Since this expression depends on St [see the right-hand side 
of equation 0 ], we simply set St to the value that yields 
d{Sr — ST^)/dt = 0. We have checked that using this ap¬ 
proximation down to z = 900 affects the power spectra by 
a fraction of a percent at most. 

For the concord ance set of cosmological parameters 
iSpergel et alJ l2003f l . with a scale-invariant primordial 
power spectrum normalized to as = 0.9 at z = 0, Fig¬ 
ure 0 compares the magnitude of the fluctuations in the 
CDM and baryon densities, and in the baryon and photon 
temperatures. For each quantity, the plot shows the dimen¬ 
sionless combination [k 3 P(k)/(2 -k 2 )] 1 ' 2 , where P(k) is the 
corresponding power spectrum of fluctuations. Note that re¬ 
gions where the fluctuations oscillate in sign (as a function 
of k) are difficult to show precisely in such a plot (e.g., the 
photons and baryons at k = 0.01-1 Mpc -1 at a = 1200, and 
the baryon temperature at k > 1000 Mpc^ 1 at 2 = 400). 
Note also that the photon temperature perturbations as 
shown are simply 1/4 of the photon density perturbations 
[see eq. 0 ], 

After recombination, two main forces affect the baryon 
density and temperature fluctuations, namely, the thermal- 


ization with the CMB and the gravitational force that at¬ 
tracts the baryons to the dark matter potential wells. As 
shown in the figure, the density perturbations in all species 
grow together [except that 5 7 = (4/3)<5b[ on scales where 
gravity is unopposed, outside the horizon (i.e., at k < 0.01 
Mpc -1 at z ~ 1000). At 2 = 1200 the perturbations in 
the baryon-photon fluid oscillate as acoustic waves on scales 
of order the sound horizon (k ~ 0.01), while smaller-scale 
perturbations in both the photons and baryons are damped 
by photon diffusion (Silk damping) and the drag of the dif¬ 
fusing photons on the baryons. Since the initial dark matter 
density perturbations increase with k, while the photon per¬ 
turbations are damped on the smallest scales by photon free 
streaming, on sufficiently small scales the power spectra of 5b 
and St roughly assume the shape of the dark matter fluctu¬ 
ation 5dm (except for the gas-pressure cutoff at the smallest 
scales), due to the effect of gravitational attraction on 5b 
and of the resulting adiabatic expansion on St- 

This evolution involves two similar physical systems. 
In each case, a target perturbation So is driven toward one 
perturbation 5i, but this forcing is opposed by coupling to 
a second perturbation 52. The values of 5i and 52 are com¬ 
parable on large scales but 1521 |5i| on small scales. As 

long as the coupling is strong, 5o ~ 82 on large scales, while 
the effect of 5i is apparent in the form of 5o on small scales 
although the coupling maintains |5o| <C |5i|. After the cou¬ 
pling weakens, the perturbation So is free to begin rising 
toward 5i, but this rise occurs only gradually. In the first 
case, 5o = 5b is driven by gravity toward 5i = 5dm, while 
mechanical coupling to 52 = 5 7 is the opposing force. In the 
second case, 5o = St is driven by adiabatic expansion toward 
5 1 = 1 5b, with resistance provided by thermal coupling to 
82 = St ■ . The mechanical coupling ends at 2 ~ 1000 while 
the thermal coupling is over by 2 ~ 200. 

The Figure also shows that the thermal tight-coupling 
approximation is accurate at the highest redshifts shown, on 
large scales since |5t — 5t 7 | -C |5t 7 | < |5b|, while on small 
scales St and 5t 7 are not strongly coupled but each is indi¬ 
vidually very small compared to 5b . Even at somewhat lower 
redshifts, 5b -C 5dm and St 5b on sub-horizon scales. By 
2 = 200 the baryon infall into the dark matter potentials is 
well advanced and adiabatic expansion is becoming increas¬ 
ingly important in setting the baryon temperature. By this 
redshift, the photon perturbations are alr eady negligible at 
k > 0.0 1 Mpc -1 , justifying their neglect bv lBarkana fc Loebl 
C3Q05c) on these scales. 


3 GROWTH OF SMALL SCALE DENSITY 
PERTURBATIONS 

On small scales (i.e., at large wavenumbers) the baryon 
perturbation growth is affected by the pressure of the gas, 
which affects the dark matter as well since the baryons con¬ 
tribute a small but significant fraction of the total gravita¬ 
tional force. The evolution of sub-horizon linear perturba¬ 
tions is described by two coupled second-order differential 
equations. The dark matter feels the combined gravity of 
itself and the baryons: 

5dm + 2R5dm = \ h 2 o ^ (/b<S b + /dm5dm) , 


(9) 
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where f2 m is the redshift zero matter density as a fraction of 
the critical density. The baryons feel both gravity and pres¬ 
sure. Prior analyses as sumed a spatially uniform baryonic 
sound speed c 3 (t) (e.g., EH fc Bertschingeill995Sh yielding 

5 b + 2 H8 b = | H 0 2 ^i (/ b 5 b + / dm 5 dm ) - ^ c 2 5 b , (10) 

2 a* a* 

where c 2 = dp/dp was calculated from the thermal evolution 
of a uniform gas undergoing Hubble expansion: 


2 = kef f _ 1 d log T \ 
p \ 3 dloga J 


( 11 ) 


where p is the mean molecular weight. This also meant that 
the gas temperature fluctuation was assumed to be propor¬ 
tional throughout space to the density fluctuation, so that 


St 

4 


ksT / p 


- 1 


( 12 ) 


In this paper we instead use the equation of state of an ideal 
gas to derive a more general equation for the baryons, 

S b +2H5 b = |lTo 2 % (f b S b + / dm 5 dm )--^ — (5 b + Sr) -(13) 

A CL CL H 

Note that during recombination, we make a similar correc¬ 
tion of the baryonic pressure force in the equations of CMB- 
FAST. In order to solve for the density perturbations, an 
evolution equation for the fluctuations in the temperature is 
also required. On sub-horizon scales and at z < 200 we can 
neglect the perturbations in the temperature and density of 
the photons in equation 0 , yielding simply 


dSr _ 2 dS b x e (t) _ 4 T 7 
~dT ~ 3 ~dt —° T dT 


(14) 


Equations 0, 0, 0, 03 and in are a closed set of 
equations describing the evolutio n of density and t e mpera- 
ture perturbations. We note that iBharadwai fc AlJ f 2004a!') 
derived a similar equation to Eq. m but solved it only for 
the case of a density perturbation that follows the Einstein- 
de Sitter growing mode 5 b oc a and thus neglected spatial 
variations in the speed of sound. 

Figure 0 shows the power spectra at redshift 20 and 
100. As time passes, the power spectrum of the baryons ap¬ 
proaches that of the dark matter except for the pressure 
cutoff, and the baryon temperature fluctuations increase as 
well. However, even during the era of the formation of the 
first galaxies (z ~ 40-20), there is still significant memory in 
the perturbations of their earlier coupling to the CMB. This 
is highlighted in Figure [3 which shows the ratios 5 b /5 dm 
and 5r/S b . In both quantities, the strong oscillations that 
are apparent at z = 400 are slowly smoothed out toward 
lower redshifts. At the largest scales, the baryons follow the 
dark matter density, and Sr/Sh evolves from 1/3 (the value 
during tight thermal coupling to the CMB) to ~ 2/3 (from 
adiabatic expansion). On smaller scales, the two ratios start 
from values -C 1 during mechanical/thermal coupling, and 
increase towards 5b/5 dm = 1 and Sr/Sh = 2/3, respectively. 
The former ratio approaches its asymptotic value earlier, 
since the baryons decouple from the photons first mechan¬ 
ically and only later thermally. At the smallest scales (be¬ 
low the baryonic Jeans scale), the baryon fluctuation is sup¬ 
pressed at all redshifts due to gas pressure, when the fc 2 
term in equation 03 dominates. It is clear from this figure 



Figure 2. Power spectra of density and temperature fluctua¬ 
tions vs. comoving wavenumber, at redshifts 100 and 20. We 
consider fluctuations in the CDM density (short-dashed curves), 
baryon density (solid curves), and baryon temperature (long- 
dashed curves). 
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Figure 3. Perturbation ratios <5b/^dm an d St/S^ vs. comoving 
wavenumber. We consider z = 400 (dotted curves), 2 : = 100 (solid 
curves), and 2 : = 20 (dashed curves). 


that the traditional assumption of 5t /£b being independent 
of scale is inaccurate at all redshifts considered here. 

Figure EJshows a detailed comparison between the fluc¬ 
tuation growth described by our improved equation 03 and 
that given by the traditional equation m In the improved 
calculation, the ratio St /5 b shows different behaviors on the 
































Growth of Linear Perturbations before the Era of the First Galaxies 5 



k [Mpc -1 ] 



k [Mpc -1 ] 


Figure 4. Perturbation ratios vs. comoving wavenumber. The up¬ 
per panel shows the ratio St /S\> • We consider the improved calcu¬ 
lation of equation IT31 at z = 100 (solid curve) and z = 20 (dotted 
curve), compared to the traditional calculation of equation m 
at z = 100 (short-dashed curve) and z = 20 (long-dashed curve). 
The lower panel shows the ratio between the perturbations in the 
improved calculation and those in the traditional (mean c s ) cal¬ 
culation. We show the ratio of 5^ values at 2 = 100 (solid curve) 
and 2 = 20 (dotted curve). Also shown is the ratio of <5<j m values 
at 2 = 100 (short-dashed curve) and 2 = 20 (long-dashed curve). 
Note that the displayed (5dm ratios have been smoothed at k > 
3000 Mpc” 1 to suppress noise due to limited numerical precision 
in CMBFAST. 


scale of the horizon, the photon acoustic oscillations, and on 
smaller scales. Although the ratio is roughly constant over 
some ranges of scales, its value differs from that of the tra¬ 
ditional calculation. The improved calculation changes the 
temperature fluctuations on all scales by > 10% at 2 = 20, 
and by much more at 2 = 100. Note that at the lower red- 
shift, St/S b at the small-fc end is higher than the adiabatic 
value of 2/3, since when T falls significantly below T 7 the 
thermal coupling tends to heat the gas more strongly in 
regions with a higher photon density [see the <5 7 term in 
eq. JHJ]. 

We also present in the lower panel of this Figure the 
ratio between the fluctuations in the improved calculation 
and those in the traditional calculation. Although the im¬ 
proved calculation has only a small effect on the dark matter 
density (< 1%) and a similarly small effect on the baryon 
density at k < 100 Mpc -1 , on smaller scales the baryon fluc¬ 
tuations are substantially affected. The baryon fluctuations 
are changed by up to 30% at 2 = 100 and 10% at 2 = 20. 
Thus, accurate initial conditions for models and simulations 
of the formation of the first galaxies require a full calcu¬ 
lation of the evolution of baryon density and temperature 
fluctuations along with the dark matter. 


4 21cm FLUCTUATIONS AND THE SPIN 
TEMPERATURE 


Quantitative calculations of 21cm absorption begin with the 
spin temperature T s , defined through the ratio between the 
number densities of hydrogen atoms, 


n i _ 9 1 c ~r«/T s 
no go 


(15) 


where subscripts 1 and 0 correspond to the excited and 
ground state levels of the 21cm transition, ( gi/go ) = 3 
is the ratio of the spin degeneracy factors of the levels, 
and T* = 0.0682K corresponds to the energy difference 
between the levels. The 21cm spin temperature is on the 
one hand radiatively coupled to the CMB temperature, 
and on the other hand coup led to the kinetic gas te m- 
perature T through collisions 

the a bsorption of Lya photons~ IWouUnivsenr^l5a ~ lFieIdl 
il95fil). For the concord ance set of cosmological parame¬ 
ters ijSnergel et al .1200311 . the mean brightness temperature 
on the sky at redsh ift 2 (re lative to the CMB itself) is 
iMadam Meikiyin &RqedliaaZl) 


T b = 28 mK 


i + 2 \ 1/2 / t s -r 7 
10 ) V T s 


xm, 


(16) 


where xm is the mean neutral fraction of hydrogen. 

In general, fluctuations in T b can be sourced by fluctua¬ 
tions in gas density, temperature, neutral fraction, radial ve¬ 
locity gradient, and Lya flux from galaxies. In this paper we 
consider the era before the formation of a significant galaxy 
population, so that a 'hi = 1 and there is no Lya flux. The 
velocity gradient term ( Sobolevl Il96fll) is in Fourier space 
jKaiserl[f987t iBharadwai fc Aliil2004bTi Sd r v T = —g 2 Sb/H, 
where d r v r is the line-of-sight gradient of the line-of-sight 
gas velocity, and /r = cos 0 k in terms of the angle 9k of k with 
respect to the line of sight. We can therefor e write t he fluctu ¬ 
ation in the brightness temperature T b as fearkana fc Loebl 
l2005al) 


S Tb (k,t) = + (35 b (k) + PtStW, (17) 


where we have defined time-dependent coefficients 0 and 
0 t [combining the relev ant explicit terms from eq. (2) of 
iBarkana fc Loeb f 2005atl 1. Before the first galaxies, these co¬ 
efficients depended only on the cosmic mean hydrogen den¬ 
sity uh and the temperatures T 7 and T. 

Figure |5] shows the redshift evolution of the various 
mean temperatures and the coefficients 0 and 0t- As dis¬ 
cussed above, after recombination the small residual frac¬ 
tion of free electrons coupled the temperature of the baryons 
to that of the photons. Thus, after recombination the lev¬ 
els of the hyperfine transition were at thermal equilibrium 
with the CMB temperature as well as the temperature of 
the gas, i.e., T s ~ T ~ T 7 , and no signal is expected. At 
2 < 200 the gas temperature dropped below the CMB tem¬ 
perature [T 7 oc (1 + 2)], eventually dropping adiabatically 
[T oc (I + 2) 2 ]. Atomic collisions in the still-dense cosmic gas 
kept T 3 below T 7 thus allowing the gas to absorb at 21cm 
against the CMB li lLoeb fc Zaldarriagal2004lf . The detection 
of 21cm fluctuations would uncover the thermal history of 
the epoch prior to the end of cosmic reionization, and mea¬ 
sure the growth of the baryon density and temperature per¬ 
turbations. 
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1 + z 

Figure 5. The upper panel plots mean temperatures vs. redshift. 
We consider T (solid curve), T 7 (short-dashed curve), and the 
21cm spin temperature T s (long-dashed curve). The lower panel 
shows the 21cm coefficients 0 (solid curve) and /3 t (dashed curve) 
vs. redshift. 


5 THERMAL SMOOTHING OF THE 21cm 
POWER SPECTRUM 


About 1% of the 21cm photons that propagate through the 
atomic medium suffer scattering by hydrogen. Scattering at 
a given observed wavelength of 21(1 + z) cm does not occur 
sharply at redshift z, but rather is distributed over a nar¬ 
row interval about «. This is due to the finite width of the 
21cm absorption line, set by a combination of the natural 
linewidth and the thermal motions of the hydrogen atoms. 
The result ing cutoff in the power s pectrum was estimated 
crudely by iBarkana fc Loen 1 200511) , but here we calculate 
it precisely along with its angular dependence. Note that 
this cutoff arises from the use of an absorption line as a 
probe, and is a separate effect from the Jeans smoothing of 
the density discussed above. 

We can obtain the observed 21cm fluctuation by start¬ 
ing with the previous, unsmoothed value and smoothing it 
along the line of sight. The cross-section for scattering has 
a thermal width A un given by 


A I'D 


2 ksT U21 
m,H c 


(18) 


where z /21 is the frequency of the center of the line. The line 
profile is given by the normalized Voigt profile function (j>, 
which represents a convolution between the natural broad¬ 
ening and the thermal broadening, and is given by 


m 


u 

7T 3 / 2 


f°° e~ v 

J-oo dv (£ - v ) 2 + « 2 


(19) 


where £ = (y — V 2 i)/Aud and u = Aio/(47tAz'd) are the 
dimensionless Voigt parameters, and A 10 = 2.87 x 10~ 15 s -1 
is the Einstein spontaneous emission coefficient for the 21cm 


transition. The line profile smoothes the absorption along 
the line-of-sight direction (comoving coordinate z c ), while no 
such smoothing occurs in the perpendicular directions ( x c 
and y c ). It is convenient to express the result as a smoothing 
with a three-dimensional window function: 

<^(r 0 ) = J $T b (r - r 0 )W(r)d 3 r . (20) 


In calculating the smoothing we translate frequency to 
comoving coordinate z c using the homogeneous relation 
{v 2 i/v) — 1 = aHzc/c, neglecting liigher-order corrections. 
The resulting window function is 

W( r ) = -^<KZ{z b ))6d{xc)6d(vc) , ( 21 ) 

where 5d is a one-dimensional Dirac delta function and we 
have defined 



k B T 

m H 


= 0.768 





kpc, (22) 


where the second expression for Rt is valid in an Einstein de- 
Sitter universe (i.e., is accurate at 2 ~ 3-150). The resulting 
power spectrum Pr b of temperature brightness fluctuations 
is related to its unsmoothed value P^ b by 

p T b = p r b \W k \ 2 ■ (23) 

The Fourier transform Wk of the window function W is eas¬ 
ily derived to be 

\ y k — e - v ^ 4 “ flT e - i ( fc /*%) 2 ( 24 ) 


where fi = cos 9k in terms of the angle 8 k of k with re¬ 
spect to the line of sight. This simple factorization of the 
Fourier transform is expected, since the line profile was con¬ 
structed by convolving a Lorentzian and a Gaussian, and so 
the compound Fourier transform is simply the product of 
the individual transforms. 

The smoothing effect is characterized by two parame¬ 
ters, Rt and u. The parameter u measures the relative im¬ 
portance of the natural compared to the thermal line width, 
and its extremely low value, u = 3.74 x 10 -19 (T / K) -1 ^ 2 , 
indicates that the natural width is in practice negligible. 
We note that collisional broadening is expected to dominate 
over the natural broadening, but even collisions with other 
atoms, which contribute the highest collision rate, produce 
broadening that is greater than the natural width by only a 
few orders of magnitude, and this is still negligible relative 
to the thermal broadening. Thus, the only significant term 
is the Gaussian term in Wk, which cuts off the 21 cm power 
spectrum at small scales, regardless of the physical source of 
the 21cm fluctuations. The cutoff drops the power spectrum 
by a factor of e at a wavevector k = 1 /(/j,Rt), where, e.g., 
Rt = 3.1 kpc at z = 100 and 1.6 kpc at 2 = 20. Larger 
cutoff scales are possible if the intergalactic medium is later 
heated by galactic radiation; e.g., T = 10 3 I< at z = 20 
implies Rt = 17 kpc. 

Si nce the unsmoothed po wer spectrum is a polynomial 
in y 2 fearkana fc Loebll2005(J) . we find 

PT b = [nP^ ( k ) + ///> (k) + P^O (fc)] e~^ RT)2 . (25) 

At each k , the power spectrum at four different values of y 
suffices to measure separately the three power spectra P^ 4 , 
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Figure 6. Power spectra of 21cm brightness fluctuations versus 
comoving wavenumber. We show the three power spectra that 
are separately observable, 4 (upper panel), 2 (middle panel), 
and P^o (lower panel). In each case we show redshifts 200, 150, 
100, 50 (solid curves, from bottom to top), 35, 25, and 20 (dashed 
curves, from top to bottom). 

P -2 . and P 0 , and the value of Rt■ If many values of /z are 
measured over a range of different fc’s (including sufficiently 
large k values), then the thermal cutoff scale Rt can be mea¬ 
sured accurately, directly yielding a combination of T and 
H. Comparison with the three power spectra can be used to 
check for self-consistency and to measure cosmological pa¬ 
rameters, since the power spectra also depend on T [through 
P and P t in eq. nn i and on the parameters (through H and 
the growth of the per t urbatio ns). 

IBarkana fc Loebl ll2005al) showed that it would be dif¬ 
ficult to measure 21cm fluctuations on the largest scales 
(k < 0.01 Mpc -1 ), si nce they are sma ll and a ngular pro¬ 
jections do not help. IBarkana fc Loebl •: 2005(h calculated 
the 21cm fluctuations on intermediate large-scale structure 
scales, so here we focus on the predicted magnitude of small- 
scale fluctuations at k > 10 Mpc -1 . Figure|S|tracks the three 
21cm power spectra from redshift 200 to 20. Each power 
spectrum rises to a maximum fluctuation level of ~ 5 mK 
at z = 50, and subsequently drops due to the reduced col- 
lisional coupling of the 21cm transition to the gas tempera¬ 
ture. Small-scale oscillations are visible at z > 100, with a 
smoother cutoff at lower redshifts. The additional thermal 
smoothing (not shown in the figure) occurs at a scale where 
the unsmoothed signal is ~ 0.1-1 mK. 


6 CONCLUSIONS 

We have computed the linear growth of fluctuations over 
the entire range of scales from those outside the horizon 
down to the scales that are affected by baryonic pressure. We 
have shown that the baryonic sound speed varies spatially, 


so that the temperature and density fluctuations must be 
tracked separately. At large wavenumbers (k > 100 Mpc -1 ) 
the growth of baryon density fluctuations is changed signifi¬ 
cantly by the inhomogeneous sound speed, by up to 30% at 
2 = 100 and 10% at z = 20. The effect on the dark matter 
evolution is much more modest since the baryons contribute 
only a small fraction of the total gravitational force felt by 
the dark matter. 

After cosmic recombination, the gas decouples mechani¬ 
cally from the photons, but remains thermally coupled down 
to z ~ 200. Starting from very low values on sub-horizon 
scales, the baryon density perturbations gradually approach 
those in the dark matter, and (somewhat later) the temper¬ 
ature perturbations approach the value expected for an adi¬ 
abatic gas. An aftermath of cosmic recombination is the sig¬ 
nature of the large-scale acoustic oscillations of the baryon- 
photon fluid at k ~ 0.01-0.2 Mpc -1 . This remnant is most 
apparent in the baryon density power spectrum but it dimin¬ 
ishes with time. At small scales, the cutoff due to gas pres¬ 
sure is apparent at all redshifts at k > 10 3 Mpc -1 , and can 
be measured directly from 21cm fluctuations which reach a 
magnitude higher than 1 mK in the redshift range z = 100- 
25. 

In addition we have calculated the smoothing of the 
21 cm power spectrum due to the finite width of the 21cm 
fine. This effect is dominated by the thermal motions of the 
scattering hydrogen atoms, and is always present regard¬ 
less of the source of the 21cm fluctuations. This effect pro¬ 
duces an anisotropic Gaussian cutoff of the power spectrum 
at small scales, k > 1 kpc -1 , similar to the characteristic 
wavenumber of suppression due to the baryonic pressure. 
Measuring the thermal smoothing along the fine of sight re¬ 
quires a high resolution in frequency, e.g., ~ 25 Hz at z = 50 
compared to the redshifted 21cm frequency of 28 MHz. 

We conclude that accurate initial conditions for analyt¬ 
ical models and numerical simulations of galaxy formation 
require a full calculation of the evolution of perturbations 
that includes the spatial fluctuations in the baryonic sound 
speed. At high redshift the gas temperature fluctuations are 
still recovering from their thermal coupling to the photons 
and are quite small at sub-horizon scales; thus, our improved 
calculation is particularly important for simulations of the 
formation of the first galaxies. 
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